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Abstract 

■ Exciting experiments in the field of atom and molecule optics have 

lately drawn much attention to the effects involved in the coherent diffrac- 
, tion of particle beams. We review the influence of the finite size of the 

particles and of their energy level spectrum on the diffraction pattern. In 
turn, we demonstrate how experimental diffraction measurements allow 
' ' , to determine these quantities of weakly bound molecules by considering 

£^ ■ the diffraction of 4 He2, 4 He3, and (02)2- 

2 ' 1 Introduction 

o : 

As early as 1927 the wave-like properties of matter were demonstrated when 
an electron beam was diffracted by a crystal of nickel [Q. Recently, propelled 
by the availability of custom-made nano-structured transmission gratings with 
periods as small as 100 nm, typical diffraction experiments from light optics 

4L ■ have been carried over to the domain of atom and molecule optics SCHEIE]- 

Nowadays, these experiments serve as precision measurement devices for, e.g., 
the interaction of atoms and molecules with solid surfaces |E1 Cj or as quantum 

qh mechanical mass spectrometers [HI E] • 

J* ' Owing to the great sensitivity of these experiments it was soon realized 

that a semi-classical theory using de Broglie waves and the results from light 
optics otherwise does not yield an accurate explanation of the experimental 



observations |§] . Only by a fully quantum mechanical approach, based on 
few-body scattering theory, could the diffraction pattern of a collimated beam 
of helium atoms and helium clusters be described ^J^J. Since all particles 
in such a beam share the same (average) velocity their de Broglie wave lengths 
are inversely proportional to their masses. Therefore, their diffraction angles are 
scaled by the same factor. On the one hand this mass-selective property allowed 
to ascribe certain experimental diffraction maxima to heavy clusters up to 4 Hc26 
|12j . On the other hand the diffraction of the diatomic helium cluster ("dimer") 
could be quantitatively studied and its bond length and binding energy could 
be determined jH]. 
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Moreover, atoms or molecules may coherently undergo internal transitions 
upon interaction with the diffracting object. Inelastic diffraction from a trans- 
mission grating involving transitions of metastable argon atoms |13| and nitro- 
gen dimers ^3] has already been observed. In these experiments, however, the 
grating served only to increase the intensity diffracted by the slits. In the case 
of coherent diffraction from many slits of a grating, inelastic diffraction gives 
rise to separate additional diffraction maxima. The diffraction angles at which 
these new maxima appear are related to the transition energies in a character- 
istic way. This was first discussed for the particularly interesting case of the 
three-body cluster of helium ( "trimer" ) m order to produce the predicted 
excited state, which is believed to be of Efimov-type (see, e. g., Ref. [Tfi]). Re- 
cently, the theory of inelastic diffraction was also applied to the examples of the 
van der Waals dimers H 2 D 2 and (D 2 )2 El- Their energy level spectra are well 
known 18 j. Therefore, they may serve to compare the theory to experimental 
data. 

2 Diffraction of very weakly bound molecules 

In this section we will review some aspects of the theory of molecular diffraction. 
Our approach has been described in detail before [11)1 111) . Therefore, we will 
focus on the final steps of the calculation which lead to the desired results. 
These results relate the experimentally recordable diffraction pattern to internal 
quantities of the molecular bound state: the energy level spectrum and the bond 
length. Only the experimentally most relevant case of diffraction from a nano- 
scale transmission grating will be considered. 

We choose the coordinate system such that the z axis is parallel to the bars 
of the grating (cf. figure As an idealization, we assume the grating to be 




Figure 1: The coordinate system is chosen such that the bars of the transmission 
grating are parallel to the z axis, and that the grating is periodic in the y 
direction. The period is denoted by d, and the slit width by s. 

translationally invariant along the z axis. Therefore, the z-component of the 
incident momentum is conserved, and the diffraction problem has effectively 
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been reduced to the xy plane. A typical transmission grating, as considered in 
this article, has a period of d = 100 nm and a slit width of s 70 nm. The 
de Broglie wave length of a helium dimer at a velocity of 1 km/s, for example, 
is AdB ~ 0.05 nm. Therefore, the diffraction angles are of the order of several 
mrad about the forward direction. 

The diffraction pattern of a well-collimated particle beam is, quite generally, 
determined by the transition amplitude between the incident and final momenta. 
For atoms, which are regarded as point-like particles (PP) here, of mass M 
with incident momentum P' and final momentum P = P' + AP, the transition 
amplitude is given by ^T] 



t pp (P P ■ P' P') 



-i^^r / dY exp(-iAP y Y/h) [l - r PF {P' x , P y ; Y)] 



(1) 



Here, r pp is the transmission function characterizing the grating. In the special 
case where the grating bars are assumed to be purely repulsive, equation l[T]l 
recovers the well-known Kirchhoff diffraction amplitude from classical optics 
PI chap. 8.5]. 

A corresponding result may be obtained for dimers. The masses of the 
constituents of a dimer are denoted by mi and m2, and the total mass is 
M = mi + ni2- The interaction between the constituents is accounted for 
by a two-body potential V(r), where r denotes the relative coordinate in the 
dimer. Bound dimer states supported by V(r) are denoted by their wave func- 
tions </> 7 (r) with corresponding binding energies Ey. In the following, we will 
always assume that the kinetic energy of an incident dimer exceeds the ground 
state energy Eq of the dimer by far. This requirement is easily met for the cases 
considered in this article. For example, the kinetic energy of a helium dimer 
at a velocity of 1 km/s is of the order of 40 meV whereas the binding energy 
of its single bound state is |2?o| ~ O.l^eV Hence upon diffraction from the 
grating, excitation as well as break-up of the dimer may occur. For an incident 
dimer in bound state </> 7 ' with momentum P', the transition amplitude to the 
state 7 with final momentum P = P' + AP is given by 

t(P x ,P v ,ct> 1 ;PL,Py,<f>Y) = 



p., 



'{2ir) 2 Mh 



J dY cxp(-iAP y Y/h) [8 iy - T^{P' x: P' y -Y)] , (2) 



where the corresponding dimer transmission function is given as the weighted 
product of two point-particle transmission functions (cf. figure EJ, 



T, 



*™{P' x ,P' y ;Y)= / d 3 r 0;(r)<Mr) (3) 



pp/mr , mi , m 2 \ PP /m 2 m 2 , 



mi 

IM' 1 '*/' 9 '' ' M a J" z IM' 1 '*!'*'' ~M V 

Finally, under similar assumptions the transition amplitude for trimers may 
be derived. For simplicity, we specialize to the case where all three constituents 
have equal masses m = M/3, as in the helium trimer 4 He3. It can be shown 
that equation J3J) also holds for trimers, if the dimer transmission function is 
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substituted by its trimer equivalent, 

*r (k- l* Y - r> + 5") T ' P (k ¥ - r- - h) ■ (4) 

The Jacobi coordinates (v,p), and an interpretation of the position arguments 
involving r y and p y are shown in figure 
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Figure 2: The position arguments of the functions r pp in equations © and 
can be interpreted as projected displacements of the positions of the constituents 
of the dimer, or trimer, from the center of mass position Y. 



Whereas the intensities of the diffraction maxima are given by the modulus 
square of the transition amplitude , two kinematic relations determine their 
angular positions. Firstly, energy conservation requires that 



(K 2 + Py)/™ + E Y = {Pi + Fj)/2M + E, 



(5) 



(for the case of atom diffraction the binding energies E 7 > , E 1 should be set to 
zero) . Secondly, the discrete periodicity of the grating implies the conservation 
of the lateral momentum P' y up to a reciprocal lattice vector, i.e. 



AP y =n2Trh/d, n = 0,±l,±2, . 



(6) 



For atoms this requirement can easily be derived from equation JQ) if the trans- 
mission function r pp (P' xl PL)Y) is assumed to have period d in its third argu- 
ment. An analogous derivation can be carried out for dimers and trimers using 
equation 

Introducing polar coordinates in the xy plane, a short calculation with equa- 
tions (0 and JHJ shows that the n-th order principal diffraction maximum is 
located at ITT] 



sin 9 r , 



E-yl \ 



-1/2 r 



P' 2 /2M ) 



2nh 



(7) 



where 9 n is related to the final momentum by P y = Psm# n , and the incident 
angle 6' is defined by P' y = P' sin#'. Two cases may be distinguished: for elastic 
(E-y = Eji) diffraction the square root factor in equation (JJJ becomes unity and 
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the classical wave optical result is recovered. For inelastic (E 1 ^ Ey) molecular 
diffraction the zeroth order principal maximum is shifted to a larger (smaller) 
angle in the case of excitation (de-excitation) , and the angular spacing of the 
diffraction maxima is increased (decreased). We note that the derivation of 
equation (JJJ) involved solely the two general conservation laws (0 and ©■ It 
applies, therefore, not only to dimers but also to larger clusters. 



3 The energy spectrum of a weakly bound molecule 



Molecular diffraction allows, in principal, to determine the transition energies 
between bound states of a weakly bound molecule. The key is to note that, 
according to equation Q , inelastic diffraction gives rise to additional diffraction 
maxima whose experimentally measurable angular spacing depends upon the 
transition energies E~ t — Ey. 

To demonstrate the results presented above we consider the deuterium mole- 
cule dimer in its ortho-ortho modification, (o-D2)2- At low nozzle temperatures 
both constituents may be assumed to be in their lowest rotor states (J = 0) 
|17| . A dimer of two such 0-D2 exhibits four rotational bound states, labeled by 
their end-over-end rotational quantum number I. All four bound states belong 
to the lowest vibrational mode |18| . The interaction with the grating induces 
transitions between these dimer states. In the special case of equal constituents, 
as in (0-02)2, however, the dimer transmission function J2J vanishes due to 
symmetry if the incoming (/') and outgoing (I) dimer states have opposite parity. 
This gives rise to the parity selection rule I' + 1 = even. 
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Figure 3: Calculated diffraction pattern of an (o-D2)2 beam at an angle of 
incidence of 8' = 14°. In-between the intense elastic maxima the weaker and 
more widely spaced inelastic maxima due to the 0— >2 transition are visible. 
Other inelastic transitions are not resolved. 



The population density of the i'-states in the incident beam may be estimated 
by their angular momentum degeneracy and the Boltzmann weight factors (21' + 
1) exp(— E11 /ksTb), where the translational beam temperature is assumed to 
be T b = 400mK [T7|. Therefore, initially only the I' = (E = -848/ueV) 
and I' = 1 (E\ = — 720/ieV) states are significantly populated J7j. Figure 
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shows a diffraction pattern for (o-D2)2 at an incident angle of 8' = 14°. The 
intense maxima are due to elastic diffraction whereas the weak maxima arise 
from the transition 0^2. The wider angular spacing of the latter is related to 
the transition energy E% — Eq 377/ieV. Inelastic maxima of other transitions 
are too weak to be resolved. 



4 The size of a weakly bound molecule 



In this section we will review how the relative intensities in an elastic molecular 
diffraction pattern may be used to determine the bond length of a molecule. 
Starting from the transition amplitudes ijTJl or l|2*)l it has been shown before 
[HI E] that the n-th order diffraction maximum intensity is of the general form 



oc 



(TTns c fi/d) + sinh 2 (imS/d) 
(im/d) 2 



exp [— (2nna / d) 2 ~\ 



(8) 



This expression contains a Kirchhoff slit function whose effective slit width s e s is 
smaller than the geometrical slit width s of the grating. The difference accounts 
for the finite size of the molecule and for its van der Waals interaction with the 
grating. (The parameters S and a in equation JHJ are also related to the van der 
Waals interaction, to irregularities of the grating, and to dimer break-up, but 
they will not be relevant in the following.) Generally the effective slit width 
assumes the velocity dependent form 



Re 



f_ s " dY [1 - r{P' x ,P' y , Y)] , v> = yJpn + Pf/M, (9) 



where r should be replaced by the respective transmission function (cf. equations 
0| and HI and Ref. for atoms). Therefore, in the case of dimers, or trimers, 
the full molecular bound state wave function enters into equation @. The 
integration extends over one slit. 

For the helium dimer, however, it was argued in Ref. [5] that the dependence 
on the full wave function may be neglected. Denoting the expectation value of 



the bond length by (r) 
approximate expression 



J d 3 r |0(r)| r one obtains, after some algebra, the 



s eff (*/)~s--y--2Re 



</2 



dY 



Y 



\{r) 



(10) 



where the momentum arguments of the transmission functions have been omit- 
ted for better readability. Figure 0] shows helium atom (circles) and dimer 
(squares) data for s c s, as determined by fitting equation (jHl to experimental 
diffraction patterns. The dashed line was calculated from equation © for he- 
lium atoms. The solid line was calculated from equation (|10|) with (r) adjusted 
such as to obtain the best agreement with the data. A short analysis reveals 
that the integrals in equations JHJ and (|10|) vanish in the limit of infinite velocity. 
Hence, the distance between the solid and the dashed lines in figure 0] allows 
to estimate the bond length to be ^(r) « 2.5 nm. A quantitative evaluation of 
the data yields the final result (r) = 5.2 ± 0.4 nm, making the helium dimer the 
largest known diatomic molecule in the ground state UJ. 
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Figure 4: Experimental values for the effective slit widths of helium atoms 
(circles) and dimers (squares). The data were obtained by fitting equation Q 
to experimental diffraction patterns measured by the Toennies group at the Max 
Planck Institut in Gottingen. The dashed line was calculated from equation @. 
The solid line was calculated from equation (|10() with (r) = 5.2 nm (best fit). 
The thin horizontal line at the top indicates the geometrical slit width s. 



For the helium trimer a similar analysis can be carried out. It follows that 
the effective slit width can be expressed in terms of the expectation value of the 
pair distance (r) — J d 3 r Jd 3 p \<fi(r, p)\ 2 r only. In this way one obtains the 
approximate formula 

8eff(l/)w*-^(r> (11) 



s/2 

-2 Re / dY 



To lessen the impact of the van der Waals interaction with the grating in favor 
of a more accurate determination of the pair distance (r) the measurements 
were carried out at an incident angle of 9' = 18°. Therefore, the geometrical 
slit width s reduces to its projection perpendicular to the incident beam. The 
data are shown in figure |SJ Analogously to the helium dimer case a comparison 
between equations (|llfl and @ shows that the distance between the dashed line 
(atoms) and the solid line (trimers) is asymptotically given by |(r). A detailed 
quantitative analysis of the pair distance is yet underway. 



5 Conclusions 

A diffraction experiment setup is a sensitive measurement apparatus for weakly 
bound molecules. We have summarized in this article how such a setup may be 
employed to study two characteristic features of the bound state of a cluster: 
The energy level spectrum and the bond length. 

In the first part we demonstrated how the energy level spectrum of a two- 
body cluster may be determined in a diffraction experiment by measuring the 
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Figure 5: Experimental values for the effective slit widths of helium atoms 
(circles) and trimers (squares) at an incident angle of 6' = 18°. The solid line 
was calculated from equation using the best fit value for (r). The thin 
horizontal line at the top indicates the projected slit width perpendicular to the 
incident beam. 



angular positions of the inelastic maxima. The results were applied to the 
realistic example of the diffraction of (o-D2)2- The derivation of the formula Q 
for the inelastic diffraction angles relies solely on two general conservation laws 
and is, therefore, also applicable to larger clusters |15j . 

In the second part we showed in which way the experimental diffraction 
intensities allow to determine the bond lengths of dimers and trimers. The 
general and intuitive concept of an effective slit width was presented, and applied 
to the clusters 4 He2 and 4 He3. A comparison with experimental data from the 
group of J. P. Toennies in Gottingen was shown. Recently, in order to determine 
the bond length of the helium trimer more accurately a non-zero incident angle 
was chosen. The analysis of this experiment is currently in progress. 
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